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We investigate theoretically and experimentally the optical beam propagation in the nematic liquid 
crystal with negative dielectric anisotropy, which is aligned homeotropically in a 80/rm-thickness 
planar cell in the presence of an externally voltage. It is predicted that the nonlocal nonlinearity 
of liquid crystal undergo an oscillatory response function with a negative nonlinear refractive index 
coefficient. We found that the oscillatory nonlocal nonlinearity can support stable bright solitons, 
which are observed in experiment. We also found that abnormal modulation instability occurs with 
infinity gain coefficient at a fixed spatial frequency, which is no depend on the beam intensity. We 
observed the modulation instability in the liquid crystal at a very low intensity (0.26IT/cm 2 ), and 
the maximum gain frequency were found kept unchange when beam power changes over 2-3 orders 
of magnitude. 


I. INTRODUCTION 


Nonlocal solitons, solitons in nonlinear media with 
a nonlocal response, have been the subject of intense 
theoretical and experimental studies in many physical 
systems [H, [If , including photorefractive crystals [3j, 
matic liquid crystals U01) lead glasses @, [§j, liq 
GS El, liquid-filled photonic crystal fibers [H, 
atomic vapors [ijJ, and Bose-Einstein condensates [IE 
El . There are various types of nonlocal solitons, such as 



vortex solitons 1 
soliton clusters 21 


and 


18[, multipole solitons [19L 
24j . Nonlocal solitons have a num¬ 
ber of interesting properties, including large phase shift 
fEt , attraction between two out-of-pliase solitons [3, 03 , 
a self-induced fractional Fourier transform |26|, and sup¬ 
pression of the collapse instability [27f ■ 

In the nonlocal response in nonlinear media, the refrac¬ 
tive index change at a given point is determined not only 
by the light intensity at that point but also by the light in¬ 
tensity near that point, which can be described by a non¬ 
local response function. There are several types of nonlo¬ 
cal response in real nonlinear media, such as the zeroth- 
order modified Bessel nonlocal response in nematic liq¬ 
uid cr yst als Q, the logarithmic nonlocal response in lead 
glass [28j, and the exponential-decay type nonlocal re¬ 
sponse in aqueous solution of rhodamine B ]29l and di¬ 
luted India ink [30j. All these response functions are 
integrable functions with non-negative value. 

Nikolov et al. found that quadratic solitons are equiv¬ 
alent to nonlocal solitons on the basis of an analogy 
between parametric interaction and diffusive nonlocal¬ 
ity [31] . The nonlocal analogy was later used to suc¬ 
cessfully describe pulse compression [32| [33 , localized 
X waves (34[, and modulational instability 35j in 
materials. Nikolov et al. showed that the second har¬ 


monic (SH), like the nonlinear refractive index forming 
a waveguide, traps the fundamental wave (FW). They 
found two types of response functions, the exponential- 
decay type and the sine-oscillatory type. On the basis 
of the sine-oscillatory type response function, one can 
understand why Buryak and Kivshar found numerically 
that quadratic solitons radiate linear waves [31, i36|. In 
2012, Esbensen et al. further investigated quadratic soli¬ 
tons with a sine-oscillatory type response function and 
found a family of analytical solutions under the strongly 
nonlocal approximation 33- However, these soliton solu¬ 
tions, whether numerical or analytical, were all unstable. 
We have shown that boundary conditions can stabilize 
this type of nonlocal soliton with sine-oscillatory response 
functions [38f. In this letter, we will demonstrate that ne¬ 
matic liquid crystal with negative dielectric anisotropy is 
new media which has oscillatory response functions. 


Nematicons, spatial optical solitons in nematic liquid 
crystals (NLC), have been the subject of intense theoret¬ 
ical and experimental studies over the past two decades 
f39l ]. The pioneering work on nematicons was reported 
in 1993 by Braun et al. [40]. They investigated the 
strong self-focusing of a laser beam in NLC in various 
geometries, from which they recognized the importance 
of molecular reorientation and anchoring at the bound¬ 
aries. Subsequently, in 1998, Warenghem et al. observed 
the beam self-trapping in capillaries filled with dye-doped 
NLC 41]. In the same year, Karpierz et al. observed the 
same phenomenon in planar cells with homeotropically 
aligned NLC [42j. They lowered the required power of 
observing nematicons to milliwatt levels. In 2000, Pec- 
cianti et al. reported on nematicon formation in pla¬ 
nar cells containing a NLC aligned homogeneously in the 
presence of an externally applied voltage [43} • They ex¬ 
tended the propagation length of nematicons for millime- 
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ter levels and found the adequate model for describing 
nematicon propagation. 

The investigation on nematicon was sparse until Conti 
et al. found that the NLC with a pretilt angle induced 
by an external low-frequency electric field is a kind of 
strongly nonlocal nonlinear medium and nematicons are 
a kind of accessible solitons 0, S S- They derived a 
simplified model and linked nematicons with quadratic 
solitons. The basic properties on nematicon have been 
revealed gradually ever since. Among others, we have to 
mention the interactions between two nematicons 00 - 
0|. Recently, Piccardi et al. reported the dark nemati¬ 
con formation in planar cells filled with dye-doped NLC 
aligned homeotropicallv, which can provide an effective 
negative nonlinearity [471 ] . The negative nonlinearity is 
realized through the guest-host interaction. 

In this letter, we observed the nematicon formation in 
planar cells containing a NLC with negative dielectric 
anisotropy and positive optical anisotropy homeotrop- 
ically aligned in the presence of an externally applied 
voltage. Following the method in 0, we got a simplified 
model with a negative Kerr coefficient and an oscillatory 
periodic response function, which can support bright ne¬ 
maticons. We outlined the connection between the sim¬ 
plified model and the equations describing quadratic soli¬ 
tons mmm- 


II. PHYSICAL MODEL AND SOLITON 

We performed a series of experiments to observe the 
nematicon formation in NLC with negative dielectric 
anisotropy. The experimental setup is illustrated in Fig. 
1(a). A light beam from a Verdi laser was focused by 
a 10X microscope objective and launched into a 80-/xm- 
thick NLC cell. The configuration of the cell was shown 
in Fig. 1(b) and the cell was filled with the homeotrop- 
ically aligned KY19-008 NLC, whose ri|| = 1.726, n± = 
1.496, average elastic constant K = 1 x 10 -11 iV, opti¬ 
cal anisotropy e° p = 0.74106, and dielectric anisotropy 
e r J = —5.3. Owing to the negative dielectric anisotropy, 
the NLC molecules will try to adjust in a low-frequency 
applied electrical field in such a manner that the molecule 
axes turn perpendicular to the direction of the electric 
field 0| . A microscope and a CCD camera were used to 
collect the light scattered above the cell during propaga¬ 
tion. 

In the presence of an externally applied (low- 
frequency) electric held E r f , the evolution of the slowly 
varying envelope A of a paraxial optical beam linearly 
polarized along X (an extraordinary light) and propagat¬ 
ing along Z can be described by the system, 

BA 

2 ik—— + V XY A + &Qe° p (sin 2 9 — sin 2 9 0 )A = 0, (1) 

O Zj 

2K (qZ2 e °( e a + e.a P ^ij-) sin(20) = 0, 

( 2 ) 



FIG. 1. Sketch of the experimental setup (a) and homeotropi- 
cally aligned nematic liquid crystal cell (b) for the observation 
of nematicons. 


where eo is the vacuum permittivity, 9 is the tilt angle 
of the NLC molecules, 9q is the nadir tilt in the absence 
of light, k = kon e (9o) with ko the vacuum wavenum¬ 
ber and n e (9o) = nj_n||/(n|cos 2 0o + n']_sin 2 9 0 ) 1 / 2 ss 

( n 2 j_ + e^sin 2 ^) 1 / 2 the refractive index of the extraordi¬ 
nary light at 9 0 , V 2 xy = d\ +d%, e r J = e || — ej_(< 0), 
e° p = n| — ri 2 j_(> 0). The term d 2 z 9 in Eq.([2|) was proven 
to be negligible compared to S7 2 XY 9, therefore it can be 
removed. The homeotropical boundaries and anchoring 
at the interfaces define 9\ x =q = 9\ x= l = tt/ 2, where L 
is the cell thickness. In the absence of light, the pretilt 
angle 9 is symmetric along X about X = L/2 (the cell 
center) and depends only on X: 

2K^ + eoe r jE 2 fS in(m = 0. ( 3 ) 

As the external bias voltage increasing, the nadir tilt an¬ 
gle 9q decreases when the voltage is above the Freedericks 
threshold. 

Furthermore, analogous to the theory for nematicon 
in positive NLC0, we can set 9 = 9 + (9/9o)A>, with $ 
being the optically induced perturbation. Noting that 
9 « 6*o and d x 9 0 in the middle of the cell when the 
beam width wq is far smaller than the cell thickness, we 
can simplify Eq. m and Eq.Q into the following system, 
which describes the coupling between A and *I>: 

BA 

2ik—— + V xy A + 2nofc 2 AnA = 0, (4) 


w^XxyAn + An = n,2\A\ 2 , (5) 

where the nonlinear changes of diffraction index An = 
e° p sin(20o) < l>/2?ro, the parameter w m ( w m > 0 for |0o| < 
7 t/ 2) is the characteristic length of the nonlinear response 
function, reads: 

w m = -^{-^-} 1/2 , ( 6 ) 

E rf e 0 \e r a f \ sin(26»o)[l - 29 0 cot(20 o )] 

and 

(C) 2 0 o sin(20 o ) 

n 2 =- ~r~ -• 

4no\e r a f \E* f [l - 2 0 o cot(20 o )] 


(7) 
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The nonlinear refractive index coefficient is defined as 
suggested by Peccianti et a ].((>'. 

The nonlocal nonlinearity in negative NLC described 
by Eq. © is extraordinary than common nonlocal non¬ 
linearity. First, the nonlinear refractive index coefficient 
n 2 is negative for negative NLC, while it is positive for 
positive NLC. Secondly and more importantly, the re¬ 
sponse function described by Eq. © is oscillatory, given 
as (in 1+2D) 



R{r) = Yo{r/w m )/(Aw 2 m ), (8) 

or R{x ) = sm(x/wm)/(2w m ) (in 1+1D), where r 2 = 
x 2 + y 2 , Y 0 (.) is zero-order Bessel function of the sec¬ 
ond kind. Then the nonlinear change of diffraction index 
An is given as 


/ OO 

R(i — r')\A(r')\ 2 d D r', 

-OO 


(9) 


where r = xi + yj for 1+2D(D=2) or r = x for 
1+1D(D=1). The response functions for negative NLC 
are not integrable, and it can not be normalized using the 
condition J_ R{r)dr = 1. In local limit case (when w m 
approach zero), we have An = ri 2 |A| 2 from Eq.© and 
nonlinearity is self-defocusing. But the response function 
Eq. © can not approach a delta function mathematically. 

It is very interesting to investigate asymptotic behav¬ 
ior of these response functions for r approach zero. From 
Eq.©, we have R(r) ~ ln(r /2w m )/2mv^ n for r —» 0, and 
it approach negative infinity. Considering that U 2 < 0 so 
the effective asymptotic function is —\n(r/2w m )/2'KW^ n , 
which is identical to the asymptotic function for posi¬ 
tive NLC. For positive NLC, the response function is 
R(r) = Ko(r/w m )/(2nw^ n ), where K 0 {.) is zero-order 
modified Bessel function of the second kind 0, 0| . The 
same asymptotic functions for positive and negative NLC 
imply that under the strongly nonlocal approximation 
( w m ^ wo), the nonlinearity in negative NLC is self- 
focusing, although the nonlinear refractive index coeffi¬ 
cient ri 2 is negative. That means we can observe the 
bright soliton in negative NLC. Our numerical simula¬ 
tion and experimental observation confirmed our predic¬ 
tion as shown in Fig.2. 

Numerical simulation of beam propagation were car¬ 
ried out based on Eq. m and Eq.© with Gaussian beam 
as an incident profile. For low power input or without 
the bias voltage, no nonlinear effect occurs and beams 
undergoe linear diffraction as shown in Fig.2(a). When 
increase beam power and bias voltage, we can found sta¬ 
ble bright solitons as shown in Fig.2(b). 

In experiment, the launched power and beam width are 
fixed to Po = 4.42mkF and wo = 4 fj,m. For rr-polarized 
beam (e-light), optically induced molecular reorientation 
occurs and contributes to the nonlinear index change. So 
soliton can be observed only for x-polarized beam. Fig. 
2(c) shows the linear diffraction in absence of the bias and 
Fig. 2(d) shows the nemation soliton formation at V = 


FIG. 2. (a)-(b) Numerical simulation of beam propagation 

based on Eq. © and Eq. © for (a) linear diffraction and (b) 
soliton propagation, (c)-(d) V-polarized e-beam propagation 
for (c) without voltage bias (linear diffraction) and (d) in the 
presence of 3.4 V at 1 kHz (soliton propagation), (e)-(f) Y- 
polarized o-beam propagation (linear diffraction) for (c) with¬ 
out voltage bias and (f) without voltage bias. Beam power 
P = 4.42mW and beam width w = 4 ym in (c)-(f). 



FIG. 3. (Color online) (a) The pretilt angle do and the Kerr 
coefficient ri 2 of the NLC vs the bias voltage V. (b) The char¬ 
acteristic length w m (a solid curve) and the critical power of a 
single soliton (circles) vs the pretilt angle do . The parameters 
are for a 80 — ym -thick cell filled with the NLC (KY19-008) 
and the critical power is a numerical result. 


3.4V, which is above the Freedericks threshold. For y- 
polarized beam (o-light), the nonlinearity does not occur 
and only diffraction beam can be observed, as shown in 
fig. 2(e) and 2(f). 

The variations of characteristic length w m and Kerr 
coefficient ri 2 dependence on pretilt angle do and the bias 
voltage V are given in Fig.3. A monotonous function of 
do on E r f (or V by V = E r fL) is described by Eq. ©. 
As shown in Fig. 3(a), 9q decreases monotonously from 
7 t/2 to 0 with increasing the bias above the Freedericks 
threshold. The Freedericks threshold can be found 0,149} 

Vf r = n( f r f, ) 1/2 - ( 1 °) 

£ 0 \£a | 

Introducing the value of e r J and K for KY19-008 NLC, 
we can get Vf r « 1.45V. For E r f higher than the 
Freedericksz threshold, the approximation 


— (EIl ) 3 

2 K E r f’ 


( 11 ) 


is satisfactory, where Ef r = Vf r /L. Therefore, we can 
clearly see from Eqs.© and © that w m and 77.2 are de¬ 
termined by V or 0q for a given NLC cell configuration. 
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FIG. 4. (Color online) (a) (1+2)D response function for neg¬ 
ative (solid) and positive (dashed) NLC. (b) The modula¬ 
tion Instability gain for different beam power, inok^n^A^ = 
1,5,10. 


As shown in Fig. 3, w m and n 2 changes nonmonotonously 
with increasing 0q and V , respectively. Both of them have 
a minimum value. There we also show the relation of the 
critical power Pq on 6$ calculating numerically based on 
Eqs. (ID ID with Gaussian beam as an incident profile. 




FIG. 5. Experimental observation of MI for different power, 
P 0 = 0.06mW, 1.0 mW, lOOmW, from top to bottom, respec¬ 
tively. The brightness of photo was increased for clear by 
software. 


III. MODULATION INSTABILITY 

The oscillatory-type response functions are singularity 
in the frequency spectrum f38|. The response function 
R(r) in the Fourier domain fcj_ is 

R{k 1 _) = l/{l-w 2 m kl), (12) 

where k± is transverse spatial frequency. It is singular 
when fcj_ = 1 jw m . This singularity induce an abnormal 
MI in negative NLC. 

Using the standard process, from Eq. [4] and [5] we find 
the MI occurs in negative NLC when 

^’c < k ± < y + y V^ 1 + 16noko\n 2 \Al), (13) 

where k c = 1 / w m is the singular frequency, Aq is the am¬ 
plitude of background planar wave. The MI gain factor 
g is given as 

n n l*±l /Ol - k c k l + An 0 kl^Alkl) 

9 {k±) =^]l- —(*rrn- (14) 

Figure 4(a) shows the profile of response functions for 
both negative and positive NLC. MI gain spectrum for 
different beam power was shown in Fig.4(b). We can see 
that at the singular frequency k c , the MI gain factor is 
also singular and it value approach to infinity. When 
beam power increases, the MI range increase but the sin¬ 
gular frequency keep unchange. We can prediction that 
MI will occur at an extremely low power due to the infin¬ 
ity value of MI gain at the singular frequency k c , which 
is independent on beam power. 

We observed the abnormal MI in negative NLC using 
the same experimental configuration as Fig. 1(a). The 
lOX-objective lens is replaced by a cylindrical lens with 
focal length 25 mm to produce an ellipse beam with 


dimension of l0g,mX2.3mm. The ellipse beam was 
launched in to the NLC cell and part of the beam was 
recorded by CCD, as shown in Fig.5. We have observed 
that MI occurs at low beam power of 0.06 mW, corre¬ 
sponding to beam intensity 0.26 W/ cm 2 . In our knowl¬ 
edge, it is the lowest intensity to induce MI in experiment. 

We read the intensity profiles from these photo and 
calculate the spatial spectrum intensities at different 21 - 
distance. By remove the initial spatial spectrum at 
z = 0 [im we get the relative MI gain spectrum for several 
distance. The maximum MI gain frequencies Km are got¬ 
ten and given in Fig.6(a). One can see that the maximum 
MI gain frequencies fluctuate around Km « 0.05^m _1 
when beam power changes from 0.06 mw to 170mIF, over 
3-order of magnitude. From our theory(fc c = \/w m ), 
Km ~ 0.05means the nonlinear characteristic 
length w m ~ 20 g,m, which is near the value gotten from 
Eq. ©. We also observed the maximum MI gain fre¬ 
quency for different bias voltage, as shown in Fig.6(a). 
At z = 400/zm, Km decreases as increasing of voltage. 
So we get the characteristic length w m decreases as in¬ 
creasing of 0 O , same as shown in Fig.3(b). However, the 
value of w m gotten from Km is always 3-4 times smaller 
than that in Fig.3(b). 

IV. CONCLUSION 

In summary, we observed experimentally the nemati- 
con formation in the planar cell containing the nematic 
liquid crystal with negative dielectric anisotropy, aligned 
homeotropically in the presence of an externally applied 
voltage. We found that the nonlocal nonlinearity in this 
system undergo an oscillatory response function with a 
negative nonlinear refractive index coefficient. The sin¬ 
gularity of the oscillatory-type response functions induce 
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FIG. 6. (a)Maximum MI gain frequency varies on the beam 
power, bias voltage V = 3.4V. (b) Maximum MI gain fre¬ 
quency varies bias voltage. 

an abnormal modulation instability, in which MI gain at 
singular frequency are infinity. In experiment, we have 
observed that MI occurs at an extremely low beam inten¬ 


sity about 0.26 W/cm 2 . So we demonstrate that negative 
NLC is new media with extraordinary oscillatory-type 
nonlocal nonlinearity. 

The singularity of the oscillatory-type response func¬ 
tion comes from the fact that the function extends to 
infinity in all space. If the oscillatory-type response func¬ 
tion is limited within a finite space, the singularity will be 
suppressed [38]. If we consider the response delay time, 
the oscillatory-type response function will be localized in 
space by a finite time. In this letter, the response delay 
is neglected, which is equivalence to that the oscillatory- 
type response function radiates for an infinity time. That 
is reason that MI gain approaches to infinity value. 

This research was supported by the National Natu¬ 
ral Science Foundation of China (Grant Nos. 11174090, 
11174091, 11074080, and 11204299). 
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